A=[1 1 1 1 1 1;
    -1 1 -1 1 -1 1;
    2 4 8 16 32 64;
    -2 4 -8 16 -32 64;
    3 9 27 81 243 729;
    4 16 64 256 1024 4096;
    ];
for i=2:6
A(:,i)=A(:,i)./factorial(i);
end
format rat
C=inv(A)
M=[1 0 0 0 0 0;0 -1 0 0 0 0;-1 0 1 0 0 0;0 1 0 -1 0 0;0 0 -1 0 1 0; 0 0 0 0 -1 1];
format rat
T5=C*inv(M)
